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Abstract. The Schwarz function has played an elegant role in understanding and 
in generating new examples of exact solutions to the Laplacian growth (or "Hclc- 
Shaw") problem in the plane. The guiding principle in this connection is the fact 
that "non-physical" singularities in the "oil domain" of the Schwarz function are 
stationary, and the "physical" singularities obey simple dynamics. We give an el- 
ementary proof that the same holds in any number of dimensions for the Schwarz 
potential, introduced by D. Khavinson and H. S. Shapiro [T7] (1989). A generaliza- 
tion is also given for the so-called "elliptic growth" problem by defining a generalized 
Schwarz potential. 

New exact solutions are constructed, and we solve inverse problems of describing 
the driving singularities of a given flow. We demonstrate, by example, how C n - 
techniques can be used to locate the singularity set of the Schwarz potential. One 
of our methods is to prolong available local extension theorems by constructing 
"globalizing families". We make three conjectures in potential theory relating to 
our investigation. 



I. Introduction 

A one-parameter family of decreasing domains, {fit}, in IR n solves the Laplacian 
growth problem with sink at x G Qt if the normal velocity, v n of the boundary 
T t := dfl t is determined by a harmonic Green's function, P(x, t), of Q t as follows. 




v n \ Tt = -VP 
AP = 0, in Q t 
P\r t = 

P(x -»■ x , t) ~ -Q ■ K(x. - x ), 



where K is the fundamental solution of the Laplace equation and Q > determines 
the suction rate at x . We can also consider Q < for the case of a source x where 
injection occurs, but this problem is stable (approaching a sphere in the limit) and is 
sometimes called the "backward-time Laplacian growth" . 

This is a nonlinear moving boundary problem, ubiquitous as an ideal model (or at 
least, first approximation) of many growth processes in nature and industry. We 
stress that we are considering here the ill-posed zero surface-tension case, where the 
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interface can encounter a cusp. The zero surface-tension case has attracted wide and 
growing attention mainly for two reasons (to be brief): (i) it has direct connections 
to many other areas such as classical potential theory, integrable systems, soliton 
theory, and random matrices (ii) it admits a miraculous complete set of explicit exact 
solutions in the two-dimensional case. 



If the domains f2 t are bounded, with Q > problem (1.1 ) actually produces a shrink- 
ing boundary. We get a growth process if Q t contains infinity so P then solves an 
exterior Dirichlet problem. In such a situation, it is common to place the sink at 
infinity by prescribing asymptotics for VP so that the flux across neighborhoods of 
infinity is proportional to Q. In the two-dimensional case this can be realized in the 
laboratory using a Hele-Shaw cell. Two sheets of glass are placed close together with 
a viscous fluid ("oil") filling the void between them. A small hole is drilled in the 
center of the top sheet and an inviscid fluid ("water") is pumped in at a constant 



rate. Then problem (1.1) serves as an ideal model for the boundary of the growing 



bubble of water. The harmonic function P(x, t), in this case, corresponds to the 



pressure in the oil domain. In other physical settings modeled by (1.1), -P(x, t) can 
be a probability, a concentration, an electrostatic field, or a temperature. Because of 
the huge amount of literature, we are limited to citing an incomplete list of papers. 
For a list of over 500 references, see fl2l . 



We are particularly attracted to this problem by the lack of explicit examples in 
dimensions higher than two. The existence, uniqueness, and regularity theory are 
well-developed in arbitrary dimensions, and in the plane there is an abundance of 
explicit, exact solutions. In dimensions higher than two, the only examples are a 
shrinking sphere (in the case when the "oil domain" Q t is bounded) or the exterior 
of a homothetically growing ellipsoid (in the case Qt is unbounded). The obvious 
explanation for this deficiency of explicit examples is a lack of conformal maps in 
higher dimensions (Liouville's Theorem) since exact solutions are usually described 
in terms of a time-dependent conformal map of the domain to the disk. However, 
exact solutions can be understood using a different tool from complex analysis, the 
Schwarz function (see j6] and Section 2 below). The following theorem relates to the 
work of S. Richardson [23] and was first stated in terms of the Schwarz function by R. 
F. Millar [2T]. Also, the discussion given by S. Howison [T3] seems to have played an 
important role in popularizing the use of the Schwarz function in studies of Laplacian 
growth. 

Theorem 1.1 (Dynamics of Singularities: M 2 ). Suppose a one-parameter family of 
domains Q t has smoothly- changing analytic boundary with Schwarz function S(z,t). 
Then it is a Laplacian growth if and only if 



(1.2) 
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The Schwarz function is only guaranteed to exist in a vicinity of a given analytic 
curve, and a priori the domain of analyticity for its time-derivative is not any better. 
Thus, it is surprising that for a Laplacian growth, ^S{z,t) coincides with a function 
analytic throughout VL t except at the singularity prescribed at the "sink". In other 
words, we can extract from equation (1.2) the following elegant description of solu- 
tions to problem (1.1 ): Singularities in Q t of the Schwarz function ofdQ t do not move 
except for one simple pole stationed at the sink x which decreases in strength at the 
rate —Q. Since equation 1.2 is given in physical coordinates rather than introducing a 
uniformized "mathematical plane", S. Howison [T3] has called it an intrinsic descrip- 
tion. In recent papers, it is typical to see a combination of the Schwarz function and 
the conformal map used to derive solutions (e.g. [2]). We will review some familiar 



examples in Section 4 and understand them completely in terms of Theorem 1.1 



The Schwarz function has been partially generalized by D. Khavinson and H. S. 
Shapiro to higher dimensions by defining a "Schwarz potential", a solution of a cer- 
tain Cauchy problem for the Laplace equation [21]. In Section 2, we will review 
the definition of the Schwarz function and the Schwarz potential before proving the 



n-dimensional version of Theorem 1.1 We also give a further generalization to the 



elliptic growth problem. The rest of the paper is guided by Theorem 2.1, which 



identifies, as the main obstacle, the problem of describing (globally) the singulari- 
ties of the Schwarz potential. In Section 4 we follow the observation made by L. 
Karp that the Schwarz potential of four-dimensional, axially-symmetric surfaces can 
be calculated exactly [15]. We give some explicit examples and also describe some 
examples of elliptic growth. In Section 5, we use C n techniques to understand the 
Schwarz potential's singularity set for a nontrivial example in IR n including the im- 
portant case n = 3. In Section 6, we discuss the connection to quadrature domains 
and Richardson's Theorem. 



2. Dynamics of Singularities 



2.1. The Schwarz Potential. Suppose T is a non-singular, real-analytic curve in 
the plane. Then the Schwarz function S(z) is the function that is complex-analytic 
in a neighborhood of T and coincides with z on T (see [6 J for a full exposition). If 
T is given algebraically as the zero set of a polynomial P(x,y), we can obtain S(z) 
by making the complex-linear change of variables z = x + iy, z = x — iy, and then 
solving for z in the equation P(^±£, = 0. For instance, suppose T is the curve 
given algebraically by the solution set of the equation (x 2 + y 2 ) 2 = a 2 (x 2 + y 2 ) +4e 2 x 2 
("C. Neumann's oval"). Then changing variables we have {zz) 2 = a 2 (zz) +e 2 (z + z) 2 . 

Solving for z gives b{z) = 2( 2 2 - £ 2 ) • 

Suppose T is more generally a nonsingular, analytic hypersurface in W 1 , and consider 
the following Cauchy problem posed in the vicinity of T. The solution exists and is 
unique by the Cauchy-Kovalevskaya Theorem. 
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(2.1) 




Definition. The solution w(x) of the Cauchy problem 2.1 is called the Schwarz Po- 
tential of r. 

Example: Let L := {x G R n : ||x|| 2 = r 2 } be a sphere of radius r. When n — 2, it is 
easy to verify that w(z) = r 2 (log \z\ + 1/2 — log(r)) solves the Cauchy Problem (2.1 ), 
and in higher dimensions the Schwarz potential is iu(x) = — ( n _2)\\x\\ n - 2 2(n-2) r ' 2 - 

In R 2 , the Schwarz function can be directly recovered from the Schwarz potential. 
Consider S(z) = 2d z w = w x — iw y . The Cauchy-Riemann equations for S follow from 
harmonicity of w, and Vw = x on F implies S(z) = z on L. 

This gives a partial generalization of the Schwarz function. The reflection principle 
associated with the Schwarz function does not generalize to higher dimensions by this 
or any other means, but the Schwarz potential retains other desirable properties. In 



particular, it allows us to generalize Theorem 1.1 to higher dimensions. 



2.2. Laplacian growth and the Schwarz potential. The following theorem gen- 



eralizes Thereom LI We consider a family of domains C W 1 with analytic bound- 
ary that also has analytic time-dependence. Such regularity assumption is natural for 
us since we are in pursuit of explicit, exact solutions. However, we should mention 
that analyticity of the boundary is a necessary condition for existence of a classical 
solution, and moreover for an analytic initial boundary there exists a unique solution 
remaining analytic with analytic time-dependence for at least some interval of time 
(see [TU] and [27]). Let w(x, t) denote the Schwarz potential of the boundary T t of 
<>/. 

Theorem 2.1 (Dynamics of Singularities: M. n ). If fit and w(x,t) are as above then 



n t solves the Laplacian growth problem (1.1) if and only if 



(2.2) |«i(x,t) = -nP(x,f) 

where n is the spatial dimension. In particular, singularities of the Schwarz potential 
in the "oil domain" do not depend on time, except for one stationed at the source 
(sink) which does not move but simply changes strength. 

Remark 1: This relates the solution of a "mathematically-posed" Cauchy problem to 
that of a "physically-posed" Dirichlet problem. 

Remark 2: Considering the relationship between the Schwarz potential and Schwarz 
function, in the case of n = 2, the Theorem says that S t = §~ t (2d z w) = -Ad z P which 
is the content of Theorem 11.11 
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Remark 3: This is closely related to the celebrated Richardson's Theorem [23J. Ac- 
tually, the connection can be established through the role that the Schwarz potential 
plays in the theory of quadrature domains (see Section [5]). Here we are able to give 
a more elementary proof consisting of two applications of the chain rule. 



Proof. Assume {&t} solves the Laplacian growth problem. We will show that for each 
t, Wti'x, t) and — nP(x, t) solve the same Cauchy problem. Then by the uniqueness 
part of the Cauchy-Kovalevskaya Theorem, they are identical. 

First, we will show that u>t(x, t)| rt = 0. Consider a point x(t) which is on T t at time 
t. The chain rule gives 

(2.3) 4w(x(t), t) = Vw(x(t), t) ■ A(t) + w t (x(t),t) 



On the other hand, by the first piece of Cauchy data in (2.1), 

d 



(2.4) 



dt 



w(x(t),t) 



d 1 



ix(f)in=x(f)-A(f) 



By the second piece of Cauchy data, 

(2.5) x(t) ■ x(t) = V«)(x(t),t) • x(t) 
Combining (2.4) and (2.5) with equation ( 2.3[ ) gives 

(2.6) ^(x,t)| rt =0 

We are done if we can show that Vw t \v t = —nVp. Given some position, x, let 
T(x) assign the value of time precisely when the boundary, T t , of the growing do- 
main passes x. Then by the Cauchy data defining the Scwharz potential (2.1), 
w Xk (xi, x 2 > x ni T(xi, X2, x n )) = Xfe. Taking the partial with respect to of 
the kth equation gives w tXk T Xk + w XkXk = 1. Summing these k equations together 
gives 



(2.7) 



Vw t • VT + Aw 



n. 



where 



Since T t is the level curve T(x) = t, VT is orthogonal to T t , and VT = p ,, 2 , 
v n is the normal velocity of Y t . Recall, v n = —VP. Thus, VT = iiypjp ■ Substitution 
into equation (2.7) gives Vw t ■ VP = — n||VP|| 2 . 



By equation (2.6), w t \r t 
Vw t \ Tt = -nVP. 



0, which implies that Vu> t and VP are parallel. So, 

□ 
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2.3. A Cauchy problem connected to Elliptic Growth. A natural generaliza- 
tion of the Laplacian growth problem is to allow a non-constant "filtration coeffi- 
cient" A and "porosity" p. Then instead of Laplace's equation the pressure satisfies 
div(ApVP) = and should have a singularity at the sink of the same type as the 
fundamental solution to this elliptic equation. Moreover, the Darcy's law determining 
the boundary velocity becomes v n = —A VP. For details, see [18], [19] . [20J. Physi- 
cally, this models the problem in a non-homogeneous medium and also relates to the 
case of Hele-Shaw cells on curved surfaces (in the absence of gravity) studied in 
Ch. 7]. 



We can formulate an equation similar to equation (2.2) that relates the pressure 
function of an elliptic growth to the time-dependence of the solution to a certain 
Cauchy problem. Let q(x) be a solution of the Poisson equation, 

(2.8) div(ApVg) = np, 

where n is the spatial dimension. Recall that a solution q can be obtained by tak- 
ing the convolution of p with the fundamental solution of the homogeneous elliptic 
equation (if one exists). We associate with an elliptic growth having filtration A and 
porosity p, the solution u of the following Cauchy problem. 



(2.9) 



div (ApVu) = near T 
u\ r = q 



Vg 



We can think of u a "generalized Schwarz potential". We have the following direct 
generalization of Theorem 2.1 As in Section 2.2, assume fl t has analytic boundary 
with analytic time-dependence. 



Theorem 2.2. IfT t = dQt an d u(x,i) is the solution to 2.9 posed on Y t then VL t is 

notion P(x, t) if am 

it(x, t) = — nP(x, t) 



an elliptic growth with pressure function P(x, t) if and only if 

d_ 
dt 



(2.10) 



Proof. As in the proof of Theorem 2.1 we show that both sides of 2.10 solve the same 
Cauchy problem. 

The first part of the argument is similar in showing that. 
(2.11) ut(x,f)| rt = 

Consider a point x(i) which is on T t at time t. The chain rule gives 
d 



(2.12) 



dt 



u(x(t), t) = Vw(x(t), t) ■ x(t) + ut(x(t),t) 
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On the other hand, by the first piece of Cauchy data in (2.9) and the chain rule again, 

(2.13) ^„(x(f), t) = j t q(x(t)) = V ff (x(0) • *(f) 
By the second piece of Cauchy data, 

(2.14) Vg(x(t)) • x(t) = Vu(x(t), t) ■ x(t) 

Combining (2.13) and ( |2.14[ ) with equation (2.12) gives the equation (2.11). 
We are done if we can show that Vut|r t = —nVP. 

We again let T(x) assign the value of time when T t passes x. Then by the Cauchy data 
defining u, V«(x, T(x)) = Vg. Multiply both sides by Xp and take the divergence: 

(2.15) XpVut • VT + div(pAVw) = div(pAVg), 
which, by definition of u and q, simplifies to 

(2.16) Wu t ■ VT = n. 

As in the proof of Theorem |2.1j , VT = ^"p i where v n is the normal velocity of T t . 
Except now v n = — AVP. Thus, VT 
gives V«i ■ VP = — n||VP|| 2 . 



-VP 
A||VP|| : 



Substitution into equation (2.16) 



By equation (2.11), u t \r t 
Vu t |r t = -nVP. 



0, which implies that Vu t and VP are parallel. So, 

□ 



Let us discuss a special case of the above. Suppose that the Problem (2.8) has a 



solution q that is entire. Let a denote Xp, and suppose a is also entire. For instance, 



A 



x 2 +l 



and p 



x 



+ 1 gives a = 1 and q 



4 i 2 
X + X . 



When a = 1 as in this 



case, the "elliptic growth" is just a Laplacian growth with a variable-coefficient law 



governing the boundary velocity. The problem (2.9) defining u becomes a Cauchy 



problem for Laplace's equation with entire data. This is the realm of the Schwarz 
potential conjecture formulated by Khavinson and Shapiro: 

Conjecture 2.3 (Khavinson, Shapiro). Suppose u solves the Cauchy problem for 
Laplace's equation posed on a nonsingular analytic surface V with real-entire data. 
Then the singularity set ofu is contained in the singularity set of the Schwarz potential 
w. 



The conjecture holds in the plane and has been shown to hold "generically" in higher 
dimensions [25J. If the conjecture is true, then for the case when a — 1, the sin- 
gularities of u are controlled throughout time by those of w. Combining this with 



Theorems 2.1 and 2.2 implies that, given a solution of the Laplacian growth problem 



(1.1), the exact same evolution can be generated amid an elliptic growth law with 
1 by a pressure function having singularities at the same locations as those of w. 



The singularities may have different time-dependence and be of a different type. 
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For instance, consider the plane and the simplest Laplacian growth of suction from 
the center of a circle so that at time t, the Schwarz function is — (a constant rate of 
suction). Let us determine the pressure required to generate the same process when 
A = an d p = x 2 + 1. To solve for u, we notice that d z u is analytic and coincides 

with 2x 3 + x on the shrinking circle. Since x = z+s ^ on the boundary, we have 
d z u = (z + ^) 3 /4 + z/2 + This is even true off the boundary since both sides are 
analytic. The singular terms are 5 ~ 8 *+ 3f2 and 1 ~ 3t ^ 2 ~ f3 ■ Thus, in order to generate 
the same "movie", the pressure must have a fundamental solution type singularity 
along with a weak "multi-pole" at the origin, both diminishing at non-constant rates. 



3. Examples 



In this section we will understand some explicit solutions in terms of the Theorems 
TTTj |2H and |2~2 



3.1. Laplacian growth two dimensions. First we review some familiar examples 
in the plane, where typically a time-dependent conformal map is introduced. Instead, 



we work entirely with the Schwarz function and check that Theorem 1.1 is satisfied. 

Example 1: Consider the family of domains D with boundary given by the curves 
{z : z = aw 2 + bw, \w\ < 1} with a, b real. The Schwarz function is given by S(z) = 
—2ab/ (a — \/a 2 + Abz) +46 3 / (a — \/a 2 + Abz) 2 which has a single- valued branch in the 
interior of the curve for appropriate parameter values a and b. The only singularities 
of the Schwarz function interior to the curve are a simple pole and a pole of order two 
at the origin. Given an initial domain from this family we can choose a one-parameter 
slice of domains so that the simple pole increases (resp. decreases) while the pole of 
order two does not change. This gives an exact solution to the Laplacian growth 
problem with injection (resp. suction) taking place at the origin. In the case of 
injection, the domain approaches a circle. In the case of suction, the domain develops 
a cusp in finite time. 



Instead of just one sink or source x with rate Q, let us extend problem 1.1 by 
allowing for multiple sinks and/or sources Xi with suction/injection rates Qi. This is 
the formulation of the problem which is often made, for instance, see the excellent 



exposition [28]. The proof of Theorem 2.1 carries through without changes so that 



the time-derivative of the Schwarz potential still coincides with -n?(x, t). The only 
difference is that now there can be multiple time-dependent point-singularities inside 

£V 

Example 2: We first consider the family of curves mentioned in Section 2.1. The 
Schwarz function of the boundary is 

_ z{a 2 + 2e 2 ) + ^V4a 4 + 4a 2 e 2 + 4e 2 z 2 
S{Z) ~ 2(z* - e 2 ) 
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which has two simple poles at z = ±e each with residue (a + 2e )/2. In order to 
satisfy the conditions imposed by Theorem 2.1 we choose e = 1 constant. Then we 
choose a(i) to be decreasing (increasing) to obtain suction (injection) at two sinks 
(sources). In the case of suction, the oval forms an indentation at the top and bottom 
and becomes increasingly pinched as the boundary approaches two tangent circles 
centered at ±1, the positions of the sinks (see Figure [T]). 

For the next example, we consider the case of Problem (1.1) where the "oil domain" 
Qt is unbounded with a sink at infinity. 

Example 3: We recall the Schwarz function for an ellipse given by the solution set of 
the equation ^ + %% = 1. Changing variables we have — = 4. Solving for 

z gives S(z) = °2^2 z + fe2 2 °a2 \/ z 2 + b 2 — a 2 . S(z) has a square root branch cut along 
the segment joining the foci ±V a 2 — b 2 , but we are only interested in the exterior 
of the ellipse, where S(z) is free of singularities. This already guarantees that any 
evolution of ellipses that has analytic time dependence can be generated by preparing 
the correct asymptotic pressure conditions to match St(z,t) which is only singular at 
infinity. In other words, we can use equation 1.2 to work backwards in specifying the 
pressure condition to generate the given flow. Since there are no finite singularities, 
we only have to specify the conditions at infinity. A realistic case is if the asymptotic 
condition is steady and isotropic: S t (z — > oo,t) ~ k/z for a constant k independent 
of t. Take a homothetic growth with a(t) = ay/t and b(t) = by/i from some initial 
ellipse with semi-axes a and b. Then S(z,t) = ^zrz z + \J ' z 2 + t(b 2 — a 2 ), and 
we have S t (z,t) = k , 1 =. where k = 2ab. 
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3.2. Examples and non-examples in M 4 . Next, we consider axially-symmetric, 
four-dimensional domains. This turns out to be simpler than the more physically 
relevant M 3 , and we will see in the next subsection that it is is equivalent to certain 
cases of elliptic growth in two and three dimensions. Lavi Karp [15] has given a 
procedure, including several explicit examples, for obtaining the singularities of the 
Schwarz potential for a domain Q that is the rotation into IR 4 of a domain in IR 2 
with Schwarz function S(z). We outline here this procedure for finding the Schwarz 
potential w{xx, x 2 , x 3 , x±). Since w solves a Cauchy problem for axially symmetric 
data posed on an axially symmetric hypersurface, with, say X\ as the axis of sym- 
metry, it is a function of two variables. Write x = x\,y = \/ x\ + x\ + x 2 , and 
w(x\, x 2 , x%, Xi) = U(x,y). What makes IR 4 convenient to work with is the fact that 
V(x,y) = y ■ U(x,y) is a harmonic function of the variables x and y. Thus, finding 
U (x, y) is reduced to solving an algebraic Cauchy problem in the plane, which can be 
done in terms of the Schwarz function S(x + iy) = S(z). The steps for writing this 
solution are outlined below. 

Step 1: Write f(z) = {S(z) ■ (S(z) - 2z). 
Step 2: Find a primitive function F(z) for f(z). 

Step 3: Write V(x,y) = Re{F(z)}. Then the Schwarz potential for Q is U(x,y) = 

V(x,y)+const. 

y 

One of the examples carried through this procedure in [15] is the family of "limacons" 
from Example 1. The result is that the Schwarz potential can be expanded about the 

origin as w(x 1 ,x 2 ,x 3 ,x 4 ) = A 2 (a,b) C^A |x|~ 2 + A x (a, 6)^|x|~ 2 + A (a, 6)|x|~ 2 + 

H(x), where H(x) is harmonic and A 2 (a,b) = -6 2 a 4 /12, A x (a,b) = ba 2 (a 2 + 26 2 )/2, 
and A (a,b) = —(a 4 + 6a 2 6 2 + 26 4 )/2. We can interpret one-parameter slices of 



this family as a Laplacian growth if we further extend problem 1.1 to allow for 
"multi-poles" (see [9 J for a discussion of multi-pole sollutions in the plane). If we 
want a Laplacian growth with just a simple sink then according to the dynamics- 



of-singularities imposed by Theorem 2.1 we need to choose the time-dependence of 
a and b so that the only singularity whose coefficient changes is the fundamental- 
solution type singularity A (a, 6)|x|~ 2 . Thus, where Ci, C 2 are constants, we need to 
have: 

( o U f A 2 (a(t),b(t)) = C 2 

[6 ' \ A 1 (a(t),b(t)) = C 1 

Unfortunately, solutions a and b of this system are locally constant so that A must 
then be constant and the whole surface does not move at all. The other examples 
of axially symmetric domains considered in [15] also require introducing multi-poles 
or even a continuum of singularities, otherwise the conditions imposed by simple 
sources/sinks leads to a similarly overdetermined system. Roughly speaking, the dif- 
ficulty is that f(z) from Step 1 above generally has more singularities than S(z). Thus, 
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if a class of domains in the plane has enough parameters to control the singularities 
and obtain a Laplacian growth, then rotation into M 4 introduces more singularities 
which must be controlled with the same number of parameters. 

There are exceptions: we can describe some exact solutions involving a simple source 
and sink (no multipoles). Consider the hypersurfaces of revolution obtained by rota- 
tion from a family of curves whose Schwarz functions have two simple poles at z = ±1 
(with not necessarily equal residues). This is a two-parameter family of surfaces; as 
parameters, we can take the residues of the Schwarz function of the profile curves. 
Let fl denote the domain in the plane bounded by the profile curve. The Schwarz 
function has the form 

S(z) = —— + + c(A, B) + d(A, B)z + z 2 H(z, A, B), 

where H(z, A, B) is analytic in Q. In what follows, we will suppress the dependence 
on A and B of higher-degree coefficients. Following Steps 1 through 3 above, we have 

f(z) = l -S{z) ■ (S(z) - 2z) 
% ( A 2 B 2 d C 2 „ , A 

2 \{z — iy (z + iy z — i z + i j 

where H\(z) is analytic in Q. Then for Step 2 we need a primitive function for f(z) 
which is 

F{z) = \ ((7^1) + (iTi) + ClLog(z C2Log(z + + H2(z) 

where H2(z) is analytic in Q. 

Then for Step 3 we have V(x,y) = Re{F(z)} = {x ^ +y2 + (a!+ f y! +y2 + C x axg(z - 1) + 
C 2 aig(z + 1) + H 3 (z). 

If we can vary A and B in a way that keeps C\ and C 2 constant, then the time- 
derivative of the Schwarz function will satisfy the dynamics-of-singularities condition. 
This seems at first to be another overdetermined problem, but actually C\ and C 2 
must be equal! Otherwise, the two branch cuts of C\ arg(^ — 1) and C 2 arg(z + 
1) will not cancel eachother outside the interval [—1,1], and the Schwarz potential 
will become singular on the surface itself. This cannot happen since the surface 
has no points (in IR n ) that are characteristic for the Cauchy problem. Therefore, 
since C\ and C 2 are equal, we spend only one dimension of our parameter space 
controlling the "non-physical" segment of singularities. This leaves freedom for the 
"physical" singularities to move, at least locally, along a one-dimensional submanifold 
of parameters. Figure [2] shows the evolution of the profile curve for a typical example 
that can be obtained in this way. 
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Figure 2. A profile of an axially-symmetric solution in IR 4 with injec- 
tion at one point and suction at another. The initial curve is plotted 
in bold. 

We omit the cumbersome formulae for the time-dependence of coefficients in the alge- 
braic description of such exact solutions. The two-parameter family of hypersurfaces 
from which they are selected can be described by the solution set of: 

+ 2(a 2 -ft 2 ) ) _ ^ ' 2(a 2 -h 2 )) 2 + X 2 + X \ + X i)J 



( 4 ( fl 2 _ tff _ a 2)( x 2 + _ / / h V 2 2 2 \ 

a 2 (a 2 -/* 2 ) " [\ Xl + 2(a 2 -h 2 )) + x 2 + x z + x *J ■ 

Similarly, one can obtain examples where the Schwarz function has three or more 
simple poles. Again the suction/injection rates will have to occur in a prescribed way 
or else the time-derivative of the Schwarz potential will have singular segments which 
are difficult to interpret physically. 

Remark: The rigidity of the inter-dependence of injection/suction rates in the above 
example is made less severe by the fact that the initial and final domains only depend 
on the total quantities injected and removed at the source and sink respectively, and 
they are independent of the rates and order of work of the source and sink (see [28] : 
the proof extends word for word to higher dimensions). Thus, injection and suction 
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can happen in any manner, say one at a time, and we will lose the "movie" but retain 
the final domain. 

In the next section we will be interested in examples that correspond to axially sym- 
metric surfaces that do not intersect the axis of symmetry. For instance, to generate 
a torus, we can choose the profile curve to be a circle of radius R and center ai, 
a > R > 0. The Schwarz function is S(z) = — ai. Step 1 gives f(z) = i/2{^-. — 

ai)(^--ai-2z). Step 2 gives F(z) = 0^ + 2R 2 aLog(z - ai) + H(z), where H(z) 

is analytic. Step 3 gives V(z) = ~^ z [ y ~^ + 2,B 2 alog \z — ai\ + R(z), where R(z) is free 

of singularities. Finally, the singular part of U(x, y) is ^^ x2 ^ 2 + 2 ^ L log \z — ai\. 

This calculation for the Schwarz potential of the four-dimensional torus was carried 
out in P and discussed in connection with a classical mean-value-property for poly- 
harmonic functions. 

3.3. Examples of elliptic growth. Examples of axially-symmetric, four-dimensional 
Laplacian growth also solve certain elliptic growth problems in two and three dimen- 
sions. The two-dimensional profile solves the planar elliptic growth problem where 
the filtration coefficient A = 1 is constant, and the porosity p(x,y) = y 2 . Indeed, we 



can check that Theorem 2.2 is satisfied. The Schwarz potential U(x,y), reduced to 



two variables, satisfies the equation AU + ^ = 0. Since div (y 2 VU) = y 2 AU + 2yU, 
then U solves the Cauchy problem 



y 



div (y 2 VU) = near T 
(3.2) { U\ r = q 

VC/| r = Vg 

with q(x,y) = (x 2 + y 2 )/8 solving the Poisson equation div (y 2 Vq) = y 2 . 

The three dimensional surfaces of revolution generated by the same profile curves solve 
a three-dimensional elliptic growth if we choose A = 1 again constant and porosity 
p(x,y,z) = a/?/ 2 + z 2 . 

It is most interesting when the domain, at least initially, avoids the line {y = 0} 
where p(x, y) vanishes. Consider, for instance, a circle of radius R centered at ai. 
This corresponds to the calculation at the end of Section 3.2| for the four-dimensional 



torus. Accordingly, a shrinking circle can be generated by a simple source combined 
with a "dipole flow" positioned at the center of the circle. 

A similar calculation applies more generally when X(x,y) = y 2 ~ m , p{x,y) = y m , with 
m a positive integer, and we can consider more general domains than circles. For 
instance, a well-known classical solution of the Laplacian growth in the plane involves 
domains VL t conformally mapped from the unit disc by polynomials. Physically the 
solution has a single sink positioned at the image of the origin under the conformal 
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map. The Schwarz function of such an tt t is meromorphic except at the sink where its 
highest order pole coincides with the degree of the polynomial. So, S(z) = ^2 i=1 fi + 
H(z), where H(z) is analytic in Q t . The solution q of div (y 2 Vq) = y m is q(x,y) = 
-. To solve for U (x, y) we first notice that V(x, y) = yU (x, y) is harmonic and 



v 



m+2 



m + 3 r ^_ 11 ^ T r ^ ym + 2 



(m+2)(m+3) 

solves a Cauchy problem with data yq(x,y) = Z +3) . Thus, d z V = — 2 r m +2) ~ 

~l (m+2)(2l) m + 2 can ^ e anar ytically continued away from the boundary. As a result, 
the flow can be generated by a combination of "multipoles" positioned at the same 
point of order not exceeding k(m + 2). This resembles the result of I. Loutsenko [19] 
stating that the same evolution can be generated by multipoles of a certain order 
under an elliptic growth where p — 1 constant and A = p^, with p a positive integer. 

This fails, in an interesting way, for negative values of m. For instance, if A(x, y) = y A 
and p(x, y) = 1/y 2 , then a circle of shrinking radius R centered at ai is not generated 
by multipoles positioned at ai. Instead, the generalized Schwarz potential U(x, y) 
has singularities at the moving point iy/a 2 — r 2 . If we instead allow the center of the 
shrinking circle to move in a way that keeps \/a 2 — r 2 constant, then the evolution 
can be generated by multi-poles at this point of order up to 3 (see figure [3]). To 
reiterate, for this evolution of shrinking circles with moving center, the generalized 
(elliptic) Schwarz potential is singular at a stationary point while the analytic Schwarz 
function has a moving singularity. Such an example has been anticipated in |18j . where 
a system of nonlinear ODEs was given governing both the strength and the moving 
position of the Schwarz function's singularities under an elliptic growth. 




-2-10 1 2 



x 

Figure 3. An elliptic growth with multi-poles of order up to 3 posi- 
tioned at z = i. The Schwarz function has a moving singularity. 
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4. The Schwarz potential in C n 

The previous sections call for a deeper look into the singularities that can arise from 
Cauchy's problem for the Laplace equation. Certain techniques can only be applied if 
the problem is "complexified" . According to the algebraic form of the initial surface 
and data, we can allow each variable to assume complex values. We then consider 
the Cauchy problem in C n where the original, physical problem becomes a relatively 
small slice. We can loosely describe the advantage of a C n -viewpoint as follows: 

Consider first the wave equation in IR n . If the initial surface is non-singular and 
algebraic and the data is real-entire, then where can the solution have singularities? 
A singularity can propagate to some point if the backwards light-cone from this point 
is tangent to the initial surface. The same is true, at least heuristically, for the 
Laplace equation, except the "light cone" emanating from a point x is the isotropic 
cone := {X^Li (zi — x°) 2 = 0}, residing in C n and only touching M. n at x . Thus, the 
initial source of the singularity is located on the part of the complexified surface only 
visible if the problem is lifted to C" . 

"Leray's principle" gives the general, precise statement of the above description of 
propagation of singularities. It is only known to be rigorously true in a neighborhood 
of the initial surface. In two dimensions, where the Schwarz potential can be calcu- 
lated easily, one can check examples to see if Leray's principle gives correct global 
results (it seems to). At the same time, this gives an appealing geometric "explana- 
tion" of the source of singularities and reveals that they are the "foci" of the curve 
in the sense of Pliicker (see [Til Section 1] and the references therein). 

In arbitrary dimensions, G. Johnsson has given a global proof [H] of Leray's principle 
for quadratic surfaces. As Johnsson points out, a major step in the proof relies on 
the fact that the gradient of a quadratic polynomial is linear, so that a certain system 
of equations can be inverted easily. This becomes much more difficult for surfaces of 
higher degree, indeed, perhaps prohibitively difficult even for specific examples. 

In this section we consider a family of surfaces of degree four, the surfaces of revolu- 



tion generated by the Neumann ovals from Example 2 in Section 3T Leray's principle 
gives an appealing geometric "explanation" for the singularities of the Schwarz po- 
tential in this example, but for the rigorous proof, we apply an ad hoc combination 
of other C n techniques (actually C 2 , after taking into account axial symmetry). 

We require the following two local extension Theorems. 

Theorem 4.1 (Zerner). Let v be a holomorphic solution of the equation Lv = 
in a domain ft C C" with C 1 boundary, and assume that the coefficients of L are 
holomorphic in Q. Let zq G dQ. If dQ is non- characteristic at z with respect to L 
then v extends holomorphically into a neighborhood of z . 

In order to define non- characteristic for a real hypersurface given by the zero set of 
0, suppose the polynomial P(x, V) expresses the leading order term of L. Then T is 
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characteristic at p if P(x, V0) vanishes at x = p. For instance, if L is the Laplacian 
then the condition for {0 = 0} to be characteristic is Y17=i fix- = 0- 
In the statement of the next theorem, M is a hypersurface (of real codimension one) 
dividing a domain Q into Q + and and X is a holomorphic hypersurface (of complex 
codimension one). 

Theorem 4.2 (Ebenfelt, Khavinson, Shapiro). Assume that M is non- characteristic 
for L at po G M, and that the holomorphic hypersurface X is non-singular at po 
and meets M transversally at that point. Then any holomorphic solution v in f2_ of 
P(Z,D)v = extends holomorphically across p^. 

Theorem 4.3. Let W(x) be the Schwarz potential of the boundary T of the domain 

Q, : = {x G E n : (ELi^i) 2 - a2 Er=i^ 2 - Ax l < °}- Then W can be analytically 
continued throughout Q\B where B is the segment {x\ G [— 1,1], = for j = 
2,..,n}. 

Remark (i): In the plane, it is easily seen that W is only singular at the endpoints 
of the segment (see Example 2 in Section 3.1). In IR 4 , it is an example done by L. 



Karp [15] , who showed that the Schwarz potential has two fundamental solution type 
singularities at the endpoints along with a uniform jump in the gradient across the 
segment. 

Remark (ii): The three-dimensional consequence of this theorem is that if we take the 



surfaces of revolution generated by the Neumann ovals in Example 2 from Section 3A_ 
then the resulting evolution is a "Laplacian growth" generated by a pressure function 
having some distribution of singularities confined to the segment {x G [— 1,1], y = 
0, z — 0}. This driving mechanism is still rather obscure though, so in the next 
section we describe an approximation by finitely many simple sinks. 

Proof. We first recall that W(x) is real- analytic in a neighborhood of each nonsingular 
point of the initial surface (in IR ra ). Indeed, if the surface is nonsingular, V$|r ^ 
so that ||V$|| 2 |r 7^ so that T is everywhere non-characteristic (in M. n ) for Laplace's 
equation and the Cauchy-Kovalevskaya Theorem applies. 

Next we write W(x±,X2, ■■■,x n ) = u(x,y) where x = x\ and y = \/ x\ + x\ + ... + x 2 w 
and we recall the axially-symmetric reduction of Laplace's equation: Au+ — ^ Uy = 0. 
Since u solves a Cauchy problem for which the data and boundary are analytic, the 
problem can be lifted to C 2 . So u(x,y) can be viewed as the restriction to M 2 of the 
solution u(X, Y), valid for X and Y each taking complex values. 

We make the linear change of variables X = Y = u zw + (n ~ 2)( " z ~" m) = 0. 
Next we make another change of variables z = f(£), w = fit]), using the conformal 
map 

m = ^ - m 
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from the unit disk to the profile of Q (for appropriate value of R) which is Neumann's 
oval (see Figure lil) . 




Figure 4. The conformal map from the disc to the Neumann oval. 
This simplifies the C 2 geometry but makes the PDE more complicated. 

Write v(£,rj) = «(/(£)> /fa))- Then v t = «*(/(£)> /fa)) ' /'(£)> and the equation 
satisfied by v is jr^ffr^ + n $I 2 m (]% - 7%) = 0, or (/(£) - f{rf))v^ + (n - 
2) (f'(r})v£ — f'(£)v n ) = 0. Upon clearing denominators, the leading-order term is 

(4.1) - a^ 2 - v 2 ) m 2 - v 2 ) - v(r 2 - e)) v, n . 

After these transformations, we arrive at a Cauchy problem posed on {£77 = 1}, with 
data v = l/2/(£)/(?y), = f(rf) /'(£)> an d ^ = /(£)/' fa)- According to the form 



of the leading-order term 4.1, the characteristic points of {£77 = 1} are (±1,±1) 
(±R,±1/R), (±1/R,±R). 

The restriction of v to the non-holomorphic set 77 = £ corresponds to the original 
problem. Since W(x) was observed to be analytic near the initial surface, v(z,w) is 
analytic in a C 2 neighborhood of the circle {£77 = 1, 77 = £}, even at the characteristic 
points (±1, ±1) on the axis of symmetry. We analytically continue v from each point 
on this circle along a radial path toward the origin. Let P{6) = (e ie ,e~ t0 ). We 
consider two cases. For the first case, 9 7^ and 7^ tt, and v can be continued up to 
the origin. For the second case, when 9 = or = n, the analytic continuation stops at 
(1/R,l/R) and (—1/R,—1/R) respectively. Thus, v can be analytically continued to 
the disk minus the segment joining these two points. This transforms (by inverting 
the conformal map) to the statement we are trying to prove about W. For each case 
we construct a globalizing family in a similar manner to the proof of the Bony-Shapira 
Theorem [4]. 

CASE 1: Suppose 9^0 and 9 7^ tt so that (e l9 , e~ %e ) is not on the axis of symmetry. 
Let < s < 1 be arbitrary. We establish the continuability of v to a neighborhood 
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of the segment {tP(9),s <t<l}. Consider the path 70 := {(re 10 , ■zrw),s < r < -} 
which is on the initial surface {£77 = 1}, and passes through none of the characteristic 
points. Thus, by the Cauchy-Kovalevskaya Theorem, v is analytic in a neighborhood 
of each point on 70. Choose E\ > small enough so that v is analytic in a E\- 
neighborhood of 70. Let flo denote this tubular (C 2 ) domain of analyticity. 

For 1 > T > s, let N £2 (T) denote the ^-neighborhood of the segment {tP(9),T < 
t < 1}. Since for each 1 > t > s, the characteristic lines through tP(6) also intersect 
70, then for a small enough e 2 , any characteristic line that intersects N £2 (T) also 



intersects f2 . Let fi^ be the set ^co(f2 U N £2 (T)) \ co(f2 )J Ufi , where co(S') denotes 
the convex hull of a set S. 

Claim. For points on dQr \ dflo, the tangent plane is a supporting hyperplane for 

proof of Claim. By definition, Qt C co(fio U N S2 (T)), and these two sets share a 
boundary near points p G dQr \ dflo. The tangent plane at p G dflr \ dQo is also a 
tangent plane for dco(QoU N £2 (T)) . By convexity, it must be a supporting hyperplane 
for co(^o U N £2 (T)). It is then also a supporting hyperplane for the subset Qt- D 

Let E := {T : v can be analytically continued to Qt}- Since 1 G £7, 22 is non-empty. 
We will show that E is both open and closed relative to [s, 1] and is therefore equal 
to [s, 1]. The fact that 22 is closed follows from the fact that the domains Qt are 
continuous and nested. To see that E is open, we apply Zerner's Theorem. Suppose 
T G E, i.e., v extends to Qt- By the Claim, the tangent plane P to Qt at p G 
<9f2T \ <9^o is a supporting hyperplane. We must have that P passes through N £2 (s). 
Otherwise, P is a supporting hyperplane for both fl and N £2 (s) and, therefore, for 
any segment joining points in each of these sets (a contradiction). Since P passes 



through N £2 (s) and not f2o; it is non-characteristic. By Theorem 4.1, v extends to a 
neighborhood of p. 

CASE 2: Suppose 9 = or 9 = ir. For specificity, say 9 = 0. Then 70 := {(r, -), s < 
passes through the characteristic point (1,1). We have already observed, 
though, that v is analytic in a neighborhood of the point (1,1). If s < 1/2?, then 
7o also passes through the characteristic points (22, 1/2?), and (1/22,2?). So, we let 
s > 1/22. Then we can still choose an e± > small enough that v is analytic in a 
^-neighborhood of 70 . We use fl again to denote this domain of analyticity. We can 
proceed in the same way as in the previous case, defining N £2 {T) and Qt, except now 
the axis of symmetry z = w intersects the advancing boundary of Qt for every value 



of T. Zerner's Theorem fails at this point of intersection, but Theorem 4^ applies 
since the complex line z = w is transversal to each of the boundaries 8Qt- Thus, we 
can again prove that the set E is open and closed relative to [s, 1], but recall that we 
assumed s > 1/22. □ 
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The method of proof can clearly be applied to other examples having axial symmetry. 
In a future study, we hope to apply C n techniques to some surfaces of degree four 
that do not have axial-symmetry. For now, we state as a conjecture what we expect 
for one such example (for simplicity we formulate it in IR 3 ). 

Conjecture 4.4. Let W{x) be the Schwarz potential of the boundary F of the domain 
Q ■= y, z) G R 3 : (x 2 + y 2 + z 2 ) 2 - (a 2 x 2 + b 2 y 2 + c 2 z 2 ) < 0} ; with a>b> c> 0. 
Then W can be analytically continued throughout Q \ B where B := {z = 0, (x 2 + 



In other words, the conjecture says that the singularity set of the Schwarz potential 
of Q, a three-dimensional version of the Neumann oval, can be confined to a set in 
the xy-p\&ne bounded by a (two-dimensional) Neumann oval. 



In order to limit the number of definitions in the exposition of our main results, we 
have so far avoided explicit mention of "quadrature domains" , but it would be remiss 
not to discuss this important connection. Also, this will allow us to give a detailed 
approximate description of the second remark made after the statement of Theorem 

M\ 

First we consider the plane. A domain Q is a quadrature domain if it admits a formula 
expressing the area integral of any analytic function / belonging to, say L 1 (f2), as a 
finite sum of weighted point evaluations of the function and its derivatives, i.e. 



where Z{ are distinct points in VL and a m k are constants independent of /. 

Suppose Q is a bounded, simply-connected domain with non-singular, analytic bound- 
ary. Then the following are equivalent. Moreover, there are simple formulas relating 
the details of each. 

(i) Q is a quadrature domain. 

(ii) The exterior logarithmic potential of Q is equivalent to that which is generated 
by finitely many interior points (allowing multipoles). 

(iii) The Schwarz function of dQ is meromorphic in Q. 

(iv) The conformal map from the disk to Q is rational. 

For the equivalence of (i) and (iii), see [HI Ch. 14]. For the equivalence of (i), (ii), 
and (iv), see [2SI Ch. 3]. 

In higher dimensions, one simply replaces "analytic" with "harmonic" in the defini- 
tion of quadrature domain. In condition (ii), "logarithmic" becomes "Newtonian". In 




5. Quadrature domains 




m=l fc=0 
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higher dimensions, "multipole" refers to a finite-order partial derivative of the funda- 
mental solution to Laplace's equation. In condition (iii), "Schwarz function" becomes 
"Schwarz potential", and instead of "meromorphic" the Schwarz potential must be 
real-analytic except for finitely many "multipoles" (as described above). Then the 
equivalence of (i), (ii), and (iii) persists in higher dimensions (see [T7J Ch. 4]). Con- 
dition (iv) of course does not extend. 

If the initial domain of a Laplacian growth is a quadrature domain, then it will 
stay a quadrature domain by virtue of the equivalence of (i) and (iii) combined with 



Theorem 2.1 Moreover, according to the formulas (omitted here) relating the details 
of (i) and (iii), the consequent time-dependence of the quadrature is the content of 
Richardson's Theorem. In the plane, the quadrature domain can be reconstructed 
from its quadrature formula, and quadrature domains are dense within natural classes 
Jordon curves; the smoother the class, the stronger the topology in which they are 
dense (see [3] and the references therein). 

Theorem 5.1 (Richardson). If Q t is a Laplacian growth with m sinks located at Xi 
with rates Qi, then for any harmonic function u 



d 



P m 

I udV = - Qiu(xj) 



If the initial domain is not a quadrature domain, then the connection of Theorem 



2.1| to Richardson's Theorem requires defining quadrature domains in the wide sense, 
allowing the quadrature formula to consist of a distribution with compact support 
contained in Q (see p2] and [21] )■ For such generalized quadrature domains, a dis- 
tribution with minimal support is called a "mother body" for the domain. The 
singularity set of the Schwarz potential gives a supporting set for the "mother body" . 

Work of Gustafsson and Sakai guarantees existence of a quadrature domain in M. n 
satisfying a prescribed quadrature formula, but besides the special examples in IR 4 
the only explicit example for n > 2 is a sphere. Moveover, little qualitative infor- 
mation is known about quadrature domains in higher dimensions besides that the 
boundary is analytic. For instance, it is not even known whether quadrature domains 
are generally algebraic (in the plane, it follows from condition (iv). We make the 
following conjecture, where we mean "quadrature domain" in the classical, restricted 
sense (otherwise the statement is trivial, since any analytic, non-singular surface is a 
quadrature domain in the wide sense): 

Conjecture 5.2. In dimensions greater than two, there exist quadrature domains 
that are not algebraic. 



For the three-dimensional example from Theorem 4.3 we were able to isolate the 
singularities for the Schwarz potential to a segment inside. Thus, Q is a quadrature 
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domain in the wide sense and has a mother body supported on this segment. We ap- 
proximate the distribution using a finite number of points on this segment. Choosing 
the points Xk = — 1 + k/2, k = 0, 1, ..,4, we numerically integrate 20 harmonic basis 
functions (writing them in terms of Legendre polynomials) over Q. If we assume a 
quadrature formula involving point evalutations at the points (xk, 0, 0), then we have 
an overdetermined linear system for the coefficients (20 equations and 5 unknowns). 
We take two surfaces, and solve the least squares problem for the coefficients (using 
the same 5 points). Then the two surfaces can be approximately described as the 
boundaries of initial and final domains driven by sinks at these points, where the 
total amount removed is given by the decrease in quadrature weight. 



Figure 5. The profile of a supposed initial (a = 1) and final (a = 
2) domain. The driving mechanism to generate the smaller domain 
starting from the larger can be approximated by certain amounts of 
suction at the indicated points. 



Suppose ^initial is given by a = 2 (see statement of Theorem 4.3 ) and fifi na i is given by 
a = 1. Then of the total volume extracted, according to the approximate description 
81% is removed at the points (±1, 0, 0), 15% at the points (±1/2, 0, 0), and 4% at the 
origin (See Figure [5]). The accuracy of this description is reflected in the fact that 
the norm of the error vector for both least squares problems is on the order of 10~ 4 . 

6. Concluding remarks 

1. The equivalent definitions of quadrature domains listed in Section [5] indicate the 
possible reformulations of the Laplacian growth problem either in terms of potential 
theory or in terms of holomorphic PDEs. The potential theory approach has attracted 
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more attention and has certain advantages such as weak formulations of Laplacian 
growth. We have focused on the holomorphic PDE approach, and in Section [4] we 
gave a glimpse of its main advantage: C n techniques. 

2. The remarks at the end of Section 2.3| mention a consequence of the Schwarz 



potential conjecture regarding Laplacian growth. It would be interesting if one could 
obtain a partial result in the other direction along the lines of "Surfaces satisfying the 
SP conjecture are preserved by Laplacian growth" . This would only be interesting in 
higher dimensions, since the conjecture is already known to be true in the plane. 



3. In Section 2.3, the discussion centered around the case when a = \p = 1 is 



constant. It is natural to consider when a is a (fixed) non-constant entire function, 
and ask if the solution q to div(aVg) = 1 generalizes the data |||x|| 2 in the Schwarz 
potential conjecture. We make the following "elliptic Schwarz potential conjecture". 

Conjecture 6.1. Suppose a is entire and that u solves the Cauchy problem on a 
nonsingular analytic surface for div{ctWu) = with entire data. Then the singularity 
set of u is contained in the singularity set of v, where v solves the Cauchy problem 
with data q the solution of div(aVq) = 1 . 

One might object to generalizing unresolved conjectures. We should point out that 
the Schwarz potential conjecture is true in the plane and simple to prove, whereas we 



do not know if Conjecture 6.1 is true in the plane. One piece of evidence for the SP 
conjecture is that the Schwarz potential developes singularities at every characteristic 
point of the initial surface [TBI Proposition 11.3]. A similar proof shows that this is 
also true for v, where {(f) = 0} being characteristic for the elliptic operator means 
Va • V0 + aV0 ■ V0 = 0. 



4. At the end of Section |3.2| we have mentioned the fact that "injection is inde- 
pendent of the order of work of sources and sinks". In other words, the Laplacian 
growths driven by different sources and sinks "commute" with eachother. We can even 
consider, say hypothetically, injection at each of infinitely many interior points of a 
domain. Then we have infinitely many processes that commute with eachother. This, 
and especially its infinitesimal version which follows from the Hadamard variational 
formula, has the form of an "integrable hierarchy". To use the preferred language 
in this setting, we have a "commuting set of flows with respect to infinitely many 
generalized times" (the "times" are the amounts that have been injected into each 
of infinitely many sources). This holds in arbitrary dimensions but has recently at- 
tracted attention in two dimensions where it is directly connected to certain integrable 
hierarchies in soliton theory (see (22] and [26]). Aspects of the higher-dimensional case 
and possible connections to other integrable systems seem completely unexplored. 

5. Quadrature domains have also appeared, often only implicitly, in solutions of 
Euler's equations. Physically, this area of fluid dynamics is much different, involving 
inviscid flow with vorticity. D. Crowdy has given a survey [5] of his own work and 
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others' (mainly in the two-dimensional case) where quadrature domains have been 
applied to vortex dynamics. 

The ellipsoid is an example of a quadrature domain in the wide sense for which 
the mother body has been calculated (see [T7J Ch. 5]). The exterior gravitational 
potential of a uniform ellipsoid coincides with that of a non-uniform density supported 
on the two-dimensional "focal ellipse" of the ellipsoid. This fact was used by Dritschel 
et al [7] as a main step in developing a model for interaction of "quasi-geostrophic" 
meteorological vortices. Actually, they didn't use the exact density of the mother 
body, but only the location of its support in order to choose a small number of 
point vortices that generate a velocity field approximating that of an ellipsoid of 
uniform vorticity. Determining the strength of the approximating point vortices is 
nothing more than interpolating the quadrature formula. Our calculation at the end 
of Section [5j and similar calculations, could have promise for extending the model in 
[7J to examples of non-ellipsoidal vortices. An important missing ingredient here is a 
stability analysis, which has been carried out for ellipsoids. 



6. Our intuition for Conjecture |5.2| is based on two suspicions regarding the axially- 
symmetric case. (1) According to the singularities of the four dimensional rotation of 



a limacon considered in Section 3.2 , the quadrature formula involves point evaluations 



up to a second-order partial derivative. On the basis of L. Karp's procedure described 



in Section |3.2[ it seems that an axially-symmetric example involving only a point 
evaluation of the function and a first-order partial with respect to x will have to be 
generated by a curve whose Schwarz function has an essential singularity at the origin. 
Then, the conformal map would be transcendental. (2) In IR 3 we expect the situation 
to be at least as bad. Following [111 Ch.s 4 and 5], one can write an integral formula 
involving a Gauss hypergeometric function for the solution of a Cauchy problem for an 
n-dimensional axially-symmetric potential. The three dimensional case of the formula 
has the same form as the four- dimensional case, except the involved hypergeometric 
function is transcendental instead of rational. 
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